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Abstract 

Higgs-boson production at the Photon Collider at TESLA is studied for masses from 
180 to 350 GeV, using realistic luminosity spectra and detector simulation. Parity of the 
SM Higgs-boson can be verified from the measurement of the interference effects in the 
W + W~ decay channel and of the angular correlations in the decays of W + W~ and ZZ 
pairs. SM-like Two Higgs Doublet Model (2HDM II) and model with the generic higgs 
couplings are considered. 



1 Introduction 



A photon collider has been proposed as a natural extension of the e + e~ linear collider project 
TESLA [1 . The physics potential of a photon collider is very rich and complementary to the 
physics program of the e + e~ and hadron-hadron colliders. It is an ideal place to study the 
mechanism of the electroweak symmetry breaking (EWSB) and the properties of the Higgs- 
boson. In paper [2j we performed realistic simulation of SM Higgs-boson production at the 
Photon Collider for W + W~ and ZZ decay channels, for Higgs-boson masses above 150 GeV. 
Due to an interference with a large Standard Model background, the process 77 — > higgs — > 
W + W~ /ZZ turns out to be sensitive not only to the Higgs-boson partial width to 77, T 77 , but 
also to the phase of the 77 — ► h coupling, 77 . A precise measurement of both T 77 and </> 77 
seems to be crucial for determination of the Higgs-boson couplings, see also [SI IH El 03 U\ ■ In 
|2] we have found that this is extremely important to combine both W + W~ and ZZ channels, 
as the first one due to a large background is is very sensitive to a phase, while the second - 
to a partial width. From combined analysis of W + W~ and ZZ decay channels the 77 partial 
width can be measured with an accuracy of 3 to 8% and the phase of the amplitude with an 
accuracy between 30 and 100 mrad. 

In this paper we extend this analysis by considering in detail various extensions of the Stan- 
dard Model, in particular the SM-like Two Higgs Doublet Model 2HDM II with and without 
CP-conservation. We perform the combined analysis of W + W~ and ZZ invariant-mass distri- 
butions and extract the corresponding Higgs-boson couplings (tan/?). For 2HDM model with 
a CP violation, we estimate precision of the H — A mixing angle. 

We also consider here a potential of establishing CP-properties of Higgs-bosons in a general 
CP-violating model from a combined analysis of the angular distributions of the W + W~ and 
ZZ decay-products. In this approach a CP properties of the Higgs boson can be studied 
independently on a production mechanism. We consider a model with generic, CP-violating 
higgs couplings to vector bosons, leading to different angular distributions for a scalar- and 
pseudoscalar-type of couplings. From such measurement the CP-parity of the observed Higgs 
state can be determined. Precision in the determination of couplings can be significantly 
improved if measurement of angular distributions is combined with the invariant-mass one. 
However, in such case additional model assumptions are needed. 



2 Measurement of the h(H) — > 77 width and phase 
from the invariant- mass distributions 

The Compton back-scattering of a laser light off high-energy electron beams is considered as 
a source of high energy, highly polarized photon beams According to the current design 
PQ, the energy of the laser photons is assumed to be fixed for all considered electron-beam 
energies. With 100% circular polarization of laser photons and 85% longitudinal polarization 
of the electron beam the luminosity spectra peaked at high 77 invariant masses is expected. 
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The analysis bases on the CompAZ parametrization jH] of the realistic luminosity spectra 
for a Photon Collider at TESLA 10J. We assume that the centre-of-mass energy of colliding 
electron beams ^/s~^, is always optimised for the production of a Higgs boson with a given 
mass. The event generation according to the cross-section formula for a vector-boson production 
including the higgs contribution [HIIUE] is done with PYTHIA 6.152 ^2j- The fast simulation 
program SIMDET version 3.01 ^3] is used to model the TESLA detector performance. The 
selection cuts are applied to select W + W~ — > qqqq and ZZ — > 11 qq events (I — //, e). 

All results presented in this paper were obtained for an integrated luminosity correspond- 
ing to one year of the photon collider running, as given by ^0]- The total photon-photon 
luminosity increases from about 600 fb _1 for ^fs^ e =305 GeV (optimal beam energy choice 



for M =200 GeV) to about 1000 fb 1 for ^fs^ e =500 GeV (optimal beam energy choice for 



M =350 GeV). 



2.1 Standard Model analysis 

In this section we summarise results of [2], where the feasibility of measuring Standard Model 
Higgs-boson production in W + W~ and ZZ decay channels in 77 option of TESLA has been 
studied for a Higgs-boson mass above 150 GeV. We have studied the signal, i.e. the Higgs-boson 
decays into the vector bosons, and the background from direct vector-bosons production. For 
the ZZ final-state a direct, i.e. non-resonant 77 — > ZZ process, is rare as it occurs via loop 
only. On contrary, the non- resonant W + W~ production is a tree-level process, and is expected 
to be large. Therefore, also an interference between the signal of W + W~ production via the 
Higgs resonance and the background from direct production may be large. This effect can 
be used to access an information about the phase of the higgs — > 77 coupling, 77 . For the 
Higgs-boson masses around 350 GeV we found that the amplitude phase 77 is more sensitive 
to the loop contributions of new, heavy charged particles than the T 77 itself. 

The invariant- mass resolution obtained from a full simulation of W + W~ and ZZ events 
based on the PYTHIA and SIMDET programs, has been parametrized as a function of the 77 
centre-of-mass energy, W^. The distributions of the reconstructed invariant mass for 77 — ► 
W + W~ and 77 — ► ZZ events are shown in Figs. ^ upper and lower panels, respectively. 
Results from a full event simulation based on PYTHIA and SIMDET are compared with the 
distribution obtained from the numerical convolution of the relevant cross-sections with the 
CompAZ spectra and the parametrized detector resolution. 

Based on a parametric description of the expected mass distributions, many experiments 
were simulated, each corresponding to one year of a Photon Collider running at TESLA at 
a nominal luminosity. The "theoretical" distributions were then fitted, simultaneously to the 
observed W + W~ and ZZ mass spectra, with the width T 77 and phase 77 considered as the only 
free parameters. Results of the fits performed for different Higgs-boson masses and at different 
electron-beam energies are shown in Figs. 121 They indicate that with a proper choice of the 
electron-beam energy, the 77 partial width can be measured with an accuracy of 3 to 8%, while 
the phase of the amplitude with an accuracy between 35 and 100 mrad, see Fig. El The 77 
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Figure 1: Distribution of the reconstructed invariant mass for 77 — > W + W~ events with a SM 
Higgs-boson mass of 180 GeV and an electron-beam energy of 152.5 GeV (upper plot) and for 
77 — ► ZZ events, with a SM Higgs-boson mass of 300 GeV and an electron-beam energy of 
250 GeV (lower plot). Results from the simulation based on PYTHIA and on the SIMDET 
detector-simulation (histogram) are compared with the distribution obtained by the numerical 
convolution of the cross-section formula with the CompAZ energy-spectra for photon-beams 
and parametrization of the detector resolution (solid line). The distribution expected without 
the higgs contribution is also shown (dashed line) 0. 
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measurement opens a new window to a precise determination of the Higgs-boson couplings and 
to search of a "new physics" . It turns out that the phase is constrained predominantly by the 
W + W~ invariant- mass distribution, thanks to large interference effects between Higgs-boson 
decay and nonresonant W + W~ production. However, two-photon width of the Higgs-boson 
is much better constrained by the measurement of the ZZ mass spectra, as the nonresonant 
background is much smaller here. A precise determination of both parameters is only possible 
when both measurements are combined. 

2.2 Determination of the Higgs-boson couplings 
in the SM-like Two Higgs Doublet Model 

In our previous analysis [2j we have considered possible deviations from the Standard Model 
predictions resulting from the loop contributions to the ^77 vertex of new heavy charged 
particles. However, deviations in the two-photon width and phase can also appear if the 
couplings of the Higgs-boson to the other particles are different than those predicted by the 
Standard Model. This possibility is studied in this paper, see also [Tlj . 

As the simplest model with non-standard higgs couplings we consider Standard Model-like 
Two Higgs Doublet Model (2HDM II) [TJ]. We study a simple version of the so called solution B, 
where the lightest Higgs-boson h couplings to fermions are fixed and they have the same values 
as in the Standard Model, except that for up-type fermions the sign is opposite. The couplings 
of h to the EW gauge-bosons may differ from the corresponding SM predictions. There is only 
one parameter in such model, (3, and all remaining couplings of h and all couplings of H and 
A can be expressed in terms of it, as presented in the table [TJ 1 

One expects significant deviations from the Standard Model predictions for a light Higgs- 
boson h, both for the two-photon width and phase, for tan/? < 1. As compared to the SM 
there is a change of a relative sign of the top-quark and the W contribution, therefore the two- 
photon width in 2HDM II (B) is significantly larger than in the Standard Model, where these 
two contributions partly cancel each other. For tan/? ~ 1 the two-photon width decreases, 

1 h and H couplings to the charged higgs boson if ± are calculated according to the 2HDM II potential .14 
assuming /x=0. 





h 


H 


A 


Xu 

Xd 
Xv 


-1 
+1 

cos(2/3) 


1 

tan /3 

— tan /3 
- sin(2/?) 


1 75 tan/3 
—i 75 tan (3 





Table 1: Couplings of the neutral Higgs-bosons to fermions and vector bosons, relative to the 
Standard Model couplings, for the considered solution B of the SM-like Two Higgs Doublet 
Model (2HDM II). 
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Figure 2: Statistical error in the determination of the Higgs-boson width T 77 (upper plot) and 
phase 77 (lower plot), from the combined fit to the observed W + W~ and ZZ mass spectra, 
as a function of the Higgs-boson mass M^. Results are given for four various centre-of-mass 
energies of colliding electron beams yfs^, as indicated in the plot. The yellow (tick) line shows 
the size of the deviations expected in the SM-like 2HDM II with an additional contribution 
due to the charged Higgs-boson of mass M H + = 800 GeV. The thin dashed line is included to 
guide the eye. [2] 
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due to the suppressed HMoop contribution (ghww ~ 0). Finally, for large values of tan/? 
(cos (2/3) w — 1) the two-photon width of the light Higgs-boson h tend to be very close to the 
expectations of the Standard Model. The only difference is due to the presence of the charged 
Higgs-boson in the loop. The difference resulting from the opposite, as compared to SM, sign 
of the down-type fermion contributions is very small and can be neglected. 

Results for the light Higgs-boson h with mass = 300 GeV, from the measurement 
of the two-photon width (times the vector-boson branching ratio) and phase are presented 
in Fig. El for various tan/5 values. Error contours (la) on the expected deviation from the 
Standard-Model predictions correspond to one year of photon collider running, i.e. L 77 rs 
840 fb _1 . In the combined fit to the invariant-mass distributions for W + W~ and ZZ events 
we take the advantage of the fact that the ratio of the corresponding branching ratios BR(h — > 
ZZ) / BR{h — > W + W~) is expected to be the same as in the Standard Model. Charged Higgs- 
boson mass is set to 800 GeV. Results presented in Fig. El show that the measurement of the 
two-photon width and phase for the light Higgs-boson h decaying into W + W~ and ZZ would 
allow a precise determination of the tan f3 value. The possible ambiguity in the measurement 
of the two-photon width is resolved by the phase measurement, which clearly distinguishes 
between low tan f3 and large tan (3 solutions. The statistical error on the extracted tan (3 value 
is shown in Fig. H|for different values of light Higgs-boson mass Mh- 

For all considered Higgs-boson masses the expected error in the tan/? determination is 
smallest for tan/5 close to 1. This is because the Higgs-boson coupling to the vector bosons is 
most sensitive to tan/3, for such values. The precision of tan/3 measurement, for tan/3 ~ 1, 
deteriorates with an increase of the Higgs-boson mass. It changes from about 1% for mass 
of 200 GeV to about 10% for mass of 350 GeV. For very high and very low tan/3, when the 
relative Higgs-boson coupling to vector bosons is close to ±1 (table [TJ), precise measurement of 
tan (3 is not possible. 

The measurement of the two-photon width and phase has been investigated also for the 
heavy scalar Higgs-boson H of the SM-like Two Higgs Doublet Model (solution B), with cou- 
plings as given in tabled In this case we consider only 0.2 < tan/3 < 1, since for tan/3 > 1 
both the top-quark and W contributions are strongly suppressed and the precision of the mea- 
surement deteriorates fast. For tan/3 ~ 1 both the two-photon width and phase of the heavy 
scalar Higgs-boson H are close to the expectations of the Standard Model (for a given M H ). 
For decreasing values of tan f3 the HMoop contribution decreases, while the top-quark one in- 
creases and as a result, the two-photon width decreases slightly for tan/3 ~ 0.5 and then starts 
to increase with tan/3. Finally, for tan/3 ~ 0.1 the Higgs-boson decays to cc start to domi- 
nate. The expected number of events with the W + W~ and ZZ decays drops rapidly and the 
measurement becomes problematic again. 

The results of the analysis of the measurement of the two-photon width and phase, for the 
heavy scalar H with mass Mh = 300 GeV, are presented in Fig. 03 Light Higgs-boson mass is 
set to 120 GeV and that of the charged Higgs-boson to 800 GeV. 

Error contours (la) on the measured deviation from the Standard Model predictions are 
presented. Results show that a precise determination of tan/3, in the considered region of 
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Figure 3: The deviation from the SM of for the light Higgs-boson h with mass Mh = 300 GeV 
in the SM-like 2HDM II (B), with charged Higgs-boson mass of 800 GeV for different values 
of tan/3. Error contours (la) on the measured deviation from the Standard Model predictions 
for the h — > 77 phase 77 and for the decay width T 77 times vector boson branching ratio 
BR{h — > VV), correspond to L 77 « 840 fb _1 . Red contour labelled 'SM' indicates the expected 
precision for the Standard Model. 

parameters, should also be possible for heavy scalar H decaying to W + W~ and ZZ. 

Shown in Fig. El is the statistical error on the extracted tan (3 values. Results are given for 
four values of heavy scalar Higgs-boson mass Mg, from 200 to 350 GeV. The expected error 
in the tan/3 determination is small (1-3 %) for low tan/3, tan/3 ^0.25. For larger values of 
tan j3, the precision depends strongly on the Higgs-boson mass. For mass of 200 GeV it changes 
between 3 and 8 %, whereas for mass of 350 GeV it is between 3 and 20 %. 



2.3 Determination of CP properties of the Higgs boson 
in the SM-like Two Higgs Doublet Model 

In the general Two Higgs Doublet Model jTSj, the mass eigenstates of the neutral Higgs-bosons 
hi, hi and h% do not match CP eigenstates h, H and A. We consider here a SM-like version of 
the CP-violating Two Higgs Doublet Model, with CP violation through a small mixing between 



8 




Figure 4: Statistical error in the determination of tan f3, for four values of light Higgs-boson 
mass M h . The simultaneous fit to the observed W + W~ and ZZ mass spectra, is considered for 
the SM-like 2HDM II (B), with charged Higgs-boson mass of 800 GeV. Centre-of-mass energy 
of colliding electron beams ^fs^e is optimised for each mass Mh- 

H and A states. We consider a simple version of such model, where the couplings of the lightest 
higgs mass-eigenstate h\ (with mass 120 GeV) are expected to correspond to the couplings of 
h boson, whereas couplings of h 2 and h 3 states can be described as the superposition of H and 
A couplings (see tabled]), as follows: 

X* 1 « Xx 

Xx ~ Xx -COS&HA + Xx -sm^HA (1) 
Xx ~ Xx- cos ®ha - Xx -sin^A 

where X denotes a fermion or a vector boson, X = u, d, V = W or Z. We study the feasibility 
of the mixing angle <&ha determination from the combined measurement of the two-photon 
width and phase for the higgs mass-eigenstate h 2 . We consider only small CP- violation effects, 
i.e. \$ HA \ < 1. 

Results of the combined analysis (WW/ZZ decay channels) of the measurement of the two- 
photon width (times vector-boson branching ratio) and phase for the scalar Higgs-boson h 2 with 
mass M h2 = 300 GeV, are presented in Fig. [7| A light Higgs-boson has mass M hl = 120 GeV 
while Mh± = 800 GeV. The simultaneous fit to the observed W + W~ and ZZ mass-spectra was 
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Figure 5: As in Fig. El for the heavy Higgs-boson H with mass 300 GeV. A light Higgs-boson 
mass is assumed to be M h = 120 GeV. 



performed. Error contours (la) on the measured deviation from the Standard Model predictions 
are shown for $^ = 0, i.e. when CP is conserved, and for CP violation with <&ha = ±0.3. 
Even a small CP-violation can significantly influence the measured two-photon width and two- 
photon phase, and therefore it is possible to determine precisely both the CP-violating mixing 
angle &ha and tan /3. 

Next, we address a question: how well one can establish the fact that in a model CP 
is conserved ? The answer can be read out from Fig. |H1 where the statistical error in the 
determination of the H — A mixing angle &ha, around $>ha — value, is shown. The results 
are presented as a function of tan/3 for four values of Higgs-boson mass Mh 2 , from 200 to 
350 GeV. As above, we assume a light Higgs-boson mass M hl = M h = 120 GeV, charged 
Higgs-boson mass of 800 GeV. 

For low value of tan/3, tan/3 pa 0.25, the expected error is of the order of 20-25 mrad 
and does not depend on the heavy Higgs-boson mass. For larger values of tan /3 the precision 
worsens for small Higgs-boson masses. For tan /3 — 1 the expected error changes from 35 mrad 
for mass of 350 GeV to about 90 mrad for mass of 200 GeV. 

To summarize this part, in the 2HDM model the Higgs-boson couplings to the vector bosons 
are always as for a scalar (say in the SM) with a modification by some functions, i.e. they are 
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proportional to g^ M . This is because there is no pseudoscalar coupling to W + W~ and ZZ in 
this model (xv = The measured W + W~ and ZZ mass-spectra are sensitive to the CP- 
violating H — A mixing angle only via a Yukawa coupling, more precisely due of the top-quark 
loop contribution to the /1277 vertex. In the next section we will consider more general case, 
with a generic pseudoscalar couplings to the gauge bosons. 



3 Secondary decay analysis 

for generic model with CP violation. 

In the general case, with a generic pseudoscalar coupling of Higgs-boson to vector bosons, not 
only the invariant-mass distribution of the reconstructed W + W~ and ZZ events, as discussed 
in sec. 2.3, but also the angular distributions of W + W~ and ZZ decay products are sensitive 
to the CP properties of the Higgs-boson (T^j. We use the following form of a coupling to vector 
boson Z of a scalar (H) and a pseudoscalar (A): 

9hzz = ig-^§- (\h ■ gT) = 9 fizz , (2) 
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Figure 7: As in FigElfor the SM-like 2HDM II (B) with CP- violation for the heavy Higgs-boson 
hi with mass 300 GeV and couplings from A light Higgs-boson has mass Mf ll = 120 GeV. 

Three values of H — A mixing angle <&ha = — 0.3, 0, 0.3 are considered. 

- in MZ (\ ?W (Pl+^)p (Pl-P2)a \ , , 

9azz - W ^— \^A ■ £ 772 ' 

cos9 w \ M| J 

where pi and P2 are the 4-momenta of the vector bosons. 

The coupling with A# corresponds to the SM-like CP-even coupling, whereas the one with 
is a general CP-odd coupling for the spin-0 boson (higgs). And we use similar expressions for 
the coupling to W, with an obvious change M z to M w in the denominator of the pseudoscalar 
coupling, and a change in normalization factor as in the SM, namely Mz/ cos 9w — > Mw- We 
see, that if A# = 1, the scalar H behaves as the SM Higgs boson. 



3.1 Angular distributions for secondary decays 

Angular variables which can be used in the analysis of CP-properties are defined in Fig. El To 
test CP-properties of the Higgs-bosons one can use the distributions of the polar angles 61 and 
02 as well as the A0 distribution, where A0 is the angle between two Z- or two W- decay 
planes. To simplify the analysis, instead of two-dimensional distribution in (cosBi, COSO2) we 
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Figure 8: Verifying a CP-conservation for the SM-like 2HDM II (B). Statistical error in the 
determination of the H — A mixing angle &ha, as a function of tan/3 value, for four values 
of heavy Higgs-boson mass Mh 2 . The simultaneous fit to the observed W + W~ and ZZ mass 
spectra, with light Higgs-boson mass of 120 GeV, charged Higgs-boson mass of 800 GeV, and 
no H - A mixing (&ha = 0), Eq. □ 

consider the distribution in a new variable, defined as 

sin 2 9i ■ sin 2 9 2 , * 

^ ~ (1 + cos 2 60 • (l + cos 2 6 2 ) ' 1 ' 

The ^-distribution corresponds to the ratio of the angular distributions expected for the decay 
of a scalar and a pseudoscalar (in a limit M h >>> M z ) [T^] . 

The angular distributions in A0 and £, expected for a scalar and a pseudoscalar higgs (Eq. EJ) 
h — > ZZ — > l + l~jj are compared in Fig. EH Both distributions clearly distinguish between 
decays of scalar and pseudoscalar higgs. For the scalar higgs the distributions are almost flat 
in both A0 and (, whereas the pseudoscalar coupling (see Eq. EJ) introduces a significant A0 
and ( dependence. From the measurement of these distributions we can try to establish the 
CP properties of the Higgs boson, even without taking into account the production mechanism. 
Such measurement should be considered as the most general one, as in principle some unknown 
heavy charged particles can significantly modify the production vertex, affecting the invariant- 
mass distributions on which the measurements of the two-photon width and two-photon phase 
are based. 
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Note, that the distribution of the higgs decay-angle 0^ is expected to be flat both for 
scalar and pseudoscalar higgses. However, the corresponding angular distribution it is not flat 
for the non-resonant background. Therefore the observed 0^ distributions, where obviously a 
signal, background and possible interference between them contribute, maybe be sensitive to 
the Higgs-boson couplings. 



3.2 Generic model with CP violating couplings 



Following the analysis described in we consider a generic higgs model with the following 
tensor couplings of a Higgs boson, denoted below by H, to ZZ and W + W~\ 

9 nzz = ig-^- { X B ■ r + A, • *~- fo+^Uf^M (5) 
cos 6^ \ Mf j 

gnww = V M W (\„ • r + X A • (Vi + P2 \^ ~ (6) 



The Standard Model couplings are reproduced for X H = 1 and = 0. Below we will 




Figure 9: The definition of the polar angles 0^, ©i and 02, and the azimuthal angles 0i 
and (f>2 for the process 77 — ► h — > ZZ — > 4 /. A</> is the angle between two Z decay planes, 
A(f) = <p2 — <pi- All polar angles are calculated in the rest frame of the decaying particle. 
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Figure 10: Normalized angular distributions in A0 (upper plot) and £ (lower plot), expected 
(Eq. 0J) for scalar and pseudoscalar higgs decays H,A-+ ZZ — > for the higgs mass of 

300 GeV. 
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consider only small deviations from the Standard Model, i.e. Xh ~ 1 and |Aa| <C 1. To simplify 
the case, we assume that higgs TC couplings to fermions are the same as in the Standard Model. 



3.3 Reconstruction of the angular distributions 

Measurement of the angular distributions has been studied using samples of ZZ events gen- 
erated with PYTHIA and SIMDET programs, as described in section I2~T1 The resolutions of 
the reconstructed decay angles 6 and 0, for leptonic and hadronic Z decays, are compared in 
Fig. ^2 f° r the scalar Higgs-boson with mass of 300 GeV. 

A very good resolution is obtained for all considered distributions, nevertheless the measured 
angular distributions are strongly affected by the selection cuts used in the analysis. Because of 
the cut on the lepton and jet angles, applied to preserve a good mass resolution (see for more 
details), we observe significant loss of the selection efficiency for events with lepton or jet emitted 
along the beam direction. Also the Durham jet algorithm used in the event reconstruction 
imposes constraints on the angular separation between leptons and jets. Both effects result in a 
highly non-uniform angular acceptance. Selection efficiencies for 77 — > ZZ — ► l + l~jj events, as 
a function of the angle A0, are presented in Fig. We consider reconstructed A0 values from 
to 7T, since we are not able to distinguish between quark and antiquark jet. The efficiencies for 
ZZ events coming from decays of scalar, pseudoscalar higgs and non-resonant background are 
compared. The polar-angle distributions differ among these three classes of events, what leads 
to a significantly different acceptances. This simulation was performed for primary electron- 
beam energy of 209 GeV optimal for a Higgs-boson mass of 300 GeV. An additional cut on the 
reconstructed ZZ invariant-mass has been introduced to optimise the signal measurement from 
the angular distribution. We found that for higgs mass of 300 GeV the accepted mass range 
lies between 286 and 312 GeV. This cuts affects mainly the non-resonant background events. 

We calculate the expected angular distributions for ZZ and W + W~ events by convo lut- 
ing the corresponding cross-section formula with the CompAZ photon-energy spectra, and in 
addition with the parametrization of the invariant-mass resolution and acceptance function 
including the angular- and jet-selection cuts. Here we neglect (small) effects of the angular 
resolution. For the measurement of the event distributions in A0 and ( we introduce the ad- 
ditional cuts on the reconstructed ZZ invariant mass (for ZZ events) or on the reconstructed 
W + W~ invariant mass as well as the higgs-decay angle ®h (for W + W~ events). The cuts were 
optimised for the smallest relative error in the signal cross-section measurement. 

Expected precision in the measurements of the angle A<p and of the variable ( distributions, 
for 77 — > ZZ — > events and 77 — > W + W~ — > 4j events are illustrated in Figs. IT3l and 

fPH respectively. Calculations were performed for primary electron-beam energy of 209 GeV 
and the Higgs-boson mass of 300 GeV. The results, presented relative to the SM prediction, 
are confronted in Figs. EH (|14|) with the expectations of the generic model with Xa_ = ±0.2 and 
Xh = 1 (A A = 0.2 and X H = 1; X A = and A^ = 1.1 ), for the reconstructed ZZ (W + W~) 
events. We see, that the generic model predicts deviations both in the normalization and in the 
shape of angular distributions. Therefor we should be able to constrain Higgs-boson couplings 
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Figure 11: Resolution in the reconstructed Z-decay angles 9 and 0, for the leptonic (©;, <fii) 
and hadronic (Qj, <pj) final states. Events were simulated with the PYTHIA and SIMDET 
programs, for a primary electron-beam energy of 209 GeV and scalar Higgs-boson mass of 
300 GeV. 



o 
c 

0) 

o 

O 0.4 



I I I I l 1 4^ 



0.2 



scalar 

pseudoscalar 
background 

! E ! 1_ 



1 



3 

A(^ zz 



Figure 12: Selection efficiency as a function of the angle A<j) between two Z decay planes, 
for ZZ — > l + l~jj events coming from the scalar higgs decays, pseudoscalar higgs decays and 
non-resonant background. Events were simulated with the PYTHIA and SIMDET programs, 
for a primary electron-beam energy of 209 GeV and Higgs-boson mass of 300 GeV. Only events 
with the reconstructed higgs mass between 286 and 312 GeV are accepted. 



17 



from the shape of the distributions, event if the overall normalization related to the Higgs-boson 
production mechanizm is not known. 



3.4 Determination of CP properties of the Higgs boson 
from the measurement of angular distributions 

Each of the considered angular distributions discussed above can be fitted with a model expec- 
tations, given in terms of the parameters and Xh- If we do not use normalisation constraints 
(i.e. do not assume a given production mechanism) there are large correlations between Xa and 
Xh in fits to a single-variable distributions, i.e. for A0, ( and Oh- Therefore we first limit 
ourselves to Xh = 1 case. Statistical errors in the determination of Xa from the fits to the 
shape of various angular distributions, to be obtained after one year of photon collider running, 
are presented in Fig. separately for the WW and ZZ channels. In the figure we present also 
the results obtained from the simultaneous (combined) fit to all three considered distributions. 
It turns out that both for ZZ and for W + W~ events the measurement of the (^-distribution 
gives the best constraints on the A^ parameter, and a precision of Xa measurement from the ( 
distribution only is close to the precision obtained from the combined fit. For Mh >250 GeV 
the expected accuracy for A^ is about 0.1 both for the ZZ and for W + W~ events. 

Finally we estimate the expected statistical errors from the combined fit to angular distribu- 
tions of ZZ and W + W~ events, they are shown in Fig. By fitting all angular distributions 
simultaneously we found the error on A^ of 0.06-0.08 for > 200 GeV. 

In case of the combined fit, correlations between Xa and Xh parameters are smaller, so both 
parameters can be fitted simultaneously. Expected statistical errors on A^ and Xh are shown 
in Fig. El For low Higgs-boson masses, Mh <250 GeV, better constrains are obtained from the 
measurement of W + W~ events, whereas for >300 GeV smaller errors are obtained from 
the measurement of ZZ events. Expected error on A^ is below 0.1 for Mh <200 GeV, while the 
corresponding error on Xh changes from about 0.1 to 0.3 with increasing Higgs-boson mass. 



3.5 Determination of CP properties from the combined 

analysis of the invariant-mass and angular distributions 

Significantly better constraints on the model parameters A^ and Xh can be obtained if we 
assume that there is no "new physics" except for the considered anomalous couplings of the 
Higgs-boson to W + W~ and ZZ. We can then calculate the expected invariant-mass distribu- 
tions for ZZ and W + W~ events. In Fig. and statistical errors on the measured invariant- 
mass distributions, expected after one year of photon collider running at a nominal luminosity, 
are compared with the predictions of the generic model. We see, that large deviations from the 
Standard Model prediction are expected for A^ 7^ 0, mainly due to the interference effects. 

Expected statistical errors in the determination of A^ and Xh from the combined fit to 
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Figure 13: Expected deviations from the Standard Model predictions for ZZ — > l + l~jj events, 
for the measurement of the angle Aifizz between two Z-decay planes (upper plot) and of the 
variable Qzz calculated from the polar angles of the Z — > and Z — > jj decays (lower plot). 
Solid (dashed) red line correspond to the model with A# = 1 and A a = 0.2 (A^ = —0.2). Signal 
and background calculations are performed for primary electron-beam energy of 209 GeV and 
the Higgs-boson mass of 300 GeV. Only events with reconstructed higgs mass between 286 and 
312 GeV are accepted. 
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Figure 14: Expected deviations from the Standard Model predictions for W + W~ — > 4j events, 
for the measurement of the angle Atfiww between two W^-decay planes (upper plot) and of the 
variable (ww calculated from the polar angles of two W — > jj decays (lower plot). Solid red 
(dashed blue) line correspond to the model with A^ = 0.2 and Xh = 1 (A a = and Xh = 1.1). 
Signal and background calculations are performed for primary electron-beam energy of 209 
GeV and the Higgs-boson mass of 300 GeV. Correlated cuts on the reconstructed mass and 
higgs decay angle were imposed to improve signal to background ratio. 



20 




200 



300 



400 

M [GeV] 



< 

< 




200 



300 



400 

M [GeV] 



Figure 15: Statistical error in the determination of A^, expected after one year of photon 
collider running, as a function of the Higgs-boson mass M. Fits were performed to the shape of 
three individual angular distributions, as indicated in the plots, separately for ZZ events (upper 
plot) and W + W~ events (lower plot). Results of the simultaneous fits to all three considered 
distributions are also given. Parameter A# is fixed to the Standard Model value, A# = 1. No 
normalisation constraints are imposed. Errors resulting from the fit were calculated for A^ = 0. 
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Figure 16: Statistical error in the determination of A^, expected after one year of photon 
collider running, as a function of the Higgs-boson mass M. Results from the simultaneous fits 
to all considered angular distributions for ZZ and W + W~ events, and from the combined fit 
for both ZZ and W + W~ events are compared. Parameter Xh is fixed to the Standard Model 
value, Xh = 1- No normalisation constraints are imposed. Errors resulting from the fit were 
calculated for A^ = 0. 

angular distributions and invariant-mass distributions are compared in Fig. EDI A relative 
normalisation of ZZ and W + W~ samples is imposed in the fit, as predicted from the cross 
section calculations. However, overall normalisation, corresponding to the integrated luminosity 
of the considered data sample, is allowed to vary. For > 200 GeV the error on A^ of about 
0.01 can be achieved from the combined fit to all considered distributions. The value of A^ is 
mainly constrained by the measurement of the invariant-mass distribution for W + W~ events. 
Expected statistical error on Xh is slightly bigger, between 0.01 and 0.02, and comes mainly 
from the measurement of the invariant-mass distribution for ZZ events. 



4 Summary 

An opportunity of measuring of the Higgs-boson properties at the the Photon Collider at 
TESLA has been studied in detail for masses between 180 and 350 GeV, using realistic lu- 
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Figure 17: As in [TBI f° r ^a (upper plot) and Ah (lower plot), for determination of both 
parameters from simultaneous fit to all considered angular distributions. 
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minosity spectra and detector simulation. We found that from the combined measurement of 
the invariant-mass distributions in the ZZ and W + W~ decay-channels, the parameters of the 
SM-like Two Higgs Doublet Model can be precisely determined. For both light and heavy scalar 
Higgs boson, a statistical precision in the determination of its coupling (tan/5) is of the order 
of 10% already after one year of Photon Collider running. In case of the Two Higgs Doublet 
Model with a CP violation, the H — A mixing angle can be determined with statistical precision 
of 20 to 90 mrad (in a small- mixing approximation). 

We also considered the combined measurement of the various angular correlations in the 
W + W~ and ZZ-decay products. For a simple generic model, with SM-couplings to fermions, 
the parameters describing a CP violation in the higgs couplings to vector bosons can be deter- 
mined to about 10%. If in addition the invariant-mass distribution is used to constrain model 
parameters, 1-2 % statistical uncertainty in the determination of the Higgs-boson couplings 
can be achieved. 
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Figure 18: Predictions for the measurement of the ZZ invariant mass Mzz, for ZZ — > l + l~jj 
events. Solid red (dashed blue) line correspond to the generic model with Aa = 0.2 and A# = 1 
(\a = and \h = 1-1)- Thin green line corresponds to the Standard Model predictions (A^ = 
and A# = 1). Points with error bars indicate the statistical precision of the measurement after 
one year of photon collider running at nominal luminosity. Signal and background calculations 
are performed for primary electron-beam energy of 209 GeV optimal for the Higgs-boson mass 
of 300 GeV. 
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Figure 19: Expected deviations from the Standard Model predictions for the measurement 
of the W + W~ invariant mass M-^Wi f° r W + W~ — *> 4j events. Invariant-mass distribution 
measured after one year of photon collider running at nominal luminosity (upper plot) and the 
ratio of this distribution to the SM predictions (lower plot) are shown. Solid red (dashed blue) 
line correspond to the model with A^ = 0.2 and X H = 1 (A^ = and X H = 1.1). Thin green line 
(upper plot) corresponds to the Standard Model predictions (A^ = and Xh = 1). Signal and 
background calculations are performed for primary electron-beam energy of 209 GeV optimal 
the Higgs-boson mass of 300 GeV. 
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Figure 20: Statistical error in the determination of A^ (upper plot) and A# (lower plot), 
expected after one year of photon collider running, as a function of the Higgs-boson mass M. 
Results from the simultaneous fits to all considered angular distributions and invariant-mass 
distributions for ZZ and W + W~ events, and from the combined fit for both ZZ and W + W~ 
events are compared. Relative normalisation of ZZ and W + W~ samples is imposed in the fit. 
Errors resulting from the simultaneous fit of both parameters were calculated for A^ = and 
Xh = 1. 
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